arXiv:1502.01689v2 [nucl-ex] 25 Jun 2015 


EUROPEAN ORGANIZATION EOR NUCLEAR RESEARCH 




RLICE 


CERN-PH-EP-2015-023 
30 January 2015 


Measurement of jet suppression 


in central Pb-Pb collisions at y^NN= 2.76 TeV 


AEICE CollaboratiorQ 


Abstract 


The transverse momentum (/jt) spectrum and nuclear modification factor (Raa) of recon¬ 
structed jets in 0-10% and 10-30% central Pb-Pb collisions at -^^nn = 2.76 TeV were mea¬ 
sured. Jets were reconstructed using the anti-/:T jet algorithm with a resolution parameter of 
R = 0.2 from charged and neutral particles, utilizing the ALICE tracking detectors and Elec¬ 
tromagnetic Calorimeter (EMCal). The jet pj spectra are reported in the pseudorapidity in¬ 
terval of |T7jet| < 0.5 for 40 < px, jet < 120GeV/c in 0-10% and for 30 < px, jet < lOOGeV/c 
in 10-30% collisions. Reconstructed jets were required to contain a leading charged parti¬ 
cle with px > 5 GeV/c to suppress jets constructed from the combinatorial background in 
Pb-Pb collisions. The leading charged particle requirement applied to jet spectra both in pp 
and Pb-Pb collisions had a negligible effect on the Raa- The nuclear modification factor 
Raa was found to be 0.28 ±0.04 in 0-10% and 0.35 ±0.04 in 10-30% collisions, indepen¬ 
dent of Px, jet within the uncertainties of the measurement. The observed suppression is in 
fair agreement with expectations from two model calculations with different approaches to 
jet quenching. 
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1 Introduction 

Diserete formulations of Quantum Chromodynamies (QCD) prediet the existenee of a eross- 
over transition from normal nuelear matter to a new state of matter ealled the Quark-Gluon 
Plasma (QGP), where the partonie eonstituents, quarks and gluons, are deeonfined. The QGP 
state is expeeted to exist at energy densities above 0.5 GeV/fm^ and temperatures above 160 
MeV[[I]|, whieh ean be reaehed in eollisions of heavy-ions at ultra-relativistie energies. The ex¬ 
istenee of the QGP is supported by the observations reported by experiments at the Relativistie 
Heavy Ion Collider (RHIC) [|2]-[5|| and at the Large Hadron Collider (LHC) [l6l-fT71. 

One way to eharaeterize the properties of the QGP is to use partons from the hard seattering 
of the partonie eonstituents in the eolliding nueleons as medium probes. Hard seattering is 
expeeted to oeeur early in the eollision evolution, produeing high transverse momentum {pj) 
partons, whieh propagate through the expanding medium and eventually fragment into jets of 
hadrons. 

Due to interaetions of the high-px partons with the medium, the energy of the partons is re- 
dueed eompared to proton-proton (pp) eollisions due to medium-indueed gluon radiation and 
eollisional energy loss (jet quenehing) [TT^fT^ . The produetion eross seetion of the initial hard 
seattered partons is ealeulable using perturbative QCD (pQCD), and the eontribution from the 
non-perturbative hadronization ean be well ealibrated via jet measurements in pp eollisions. 

Jet quenehing has been observed at RHIC [l20] - l^ and at the LHC dH [IB UTl l30] - l4T]l via the 
measurement of inelusive hadron and jet produetion at high pi, di-hadron angular eorrelations 
and the dijet energy imbalanee. In all eases, the measured observable is found to be strongly 
modified in eentral heavy-ion eollisions relative to pp eollisions, when eompared to expeetations 
based on treating heavy-ion eollisions as an ineoherent superposition of independent nueleon- 
nueleon eollisions. 

Measurements of the jet kinematies are expeeted to be more elosely eorrelated to the initial par- 
ton kinematies than measurements of high-px hadrons. Jets are usually reeonstrueted by group¬ 
ing measured partieles within a given distanee, e.g. a eone with radius R. The interaetion with 
the medium ean result in a broadening of the jet shape, a softening of the jet fragmentation [l42l 
leading to an inerease of out-of-eone gluon radiation [|^ with respeet to jets reeonstrueted in 
pp eollisions liTTlI . Therefore, for a given jet resolution parameter R and a fixed initial parton 
energy, the energy of jets reeonstrueted in heavy-ion eollisions is expeeted to be smaller than 
those reeonstrueted in pp eollisions. 

Jet measurements in heavy-ion eollisions are ehallenging sinee a single event ean have multiple, 
possibly overlapping, jets from independent nueleon-nueleon seatters, as well as eombinatorie 
“jets” from the large, partially eorrelated and fluetuating baekground of low transverse mo¬ 
mentum partieles. Consequently, jet reeonstruetion in heavy-ion eollisions requires a robust 
jet-signal definition, and a proeedure to eorreet for the presenee of the large baekground and its 
assoeiated region-to-region fluetuations . 

The results reported in this letter are from lead-lead (Pb-Pb) eollision data at an energy per 
nueleon pair of = 2.76 TeV reeorded by the ALICE deteetor in 2011. Charged partieles 
are reeonstrueted with the Inner Traeking System (ITS) and the Time Projeetion Chamber (TPC) 
down to pt of 0.15 GeV/c. Neutral partieles, exeluding neutrons and ^^s, are reeonstrueted with 
the Eleetromagnetie Calorimeter (EMCal) down to a transverse energy of the EMCal elusters 
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of 0.3 GeV. For jet reeonstruetion, we followed the approaeh applied in Refs. [|451|43, where 
the average energy density of the event was subtraeted from the signal jets on a jet-by-jet basis, 
and the deteetor and baekground effeets were eorreeted on an ensemble basis via an unfolding 
proeedure. The signal jets were obtained using the anti-kx jet algorithm ll47l with a resolution 
parameter of R = 0.2 in the pseudorapidity range of |T]jet| < 0.5. Signal jets were required 
to eontain at least one eharged partiele with pr > 5 GeV/c. The eorreeted jet pj speetra and 
nuelear modifieation faetors (Raa) are reported for 40 < px, jet < 120 GeV/c in 0-10% and for 
30 < px, jet <100 GeV/c in 10-30% eentral Pb-Pb eollisions and the eorreeted jet px speetrum 
for 20 < Px, jet < 120 GeV/c in pp eollisions at a/s' = 2.76 TeV from 13.6 nb^^ reeorded in 
2011. The Raa is eompared to expeetations from two jet quenehing model ealeulations with 
different approaehes, deseribed later, in order to test the sensitivity of the observable to the 
energy density via the eentrality dependenee, and to the parton energy seale via the momentum 
dependenee. 

2 Experimental setup 

For a eomplete deseription of the ALICE deteetor and its performanee see Refs. Il48l and fl9l . 
respeetively. The analysis presented here relies mainly on the ALICE traeking system and 
EMCal, both of whieh are loeated inside a large solenoidal magnet with field strength 0.5 T. 

The traeking system eonsists of the ITS, a high-preeision six-layer silieon deteetor system with 
the inner layer at 3.9 em and the outer at 43 em from the eenter of the deteetor, and the TPC 
with a radial extent of 85-247 em, provides up to 159 independent spaee points per traek. The 
two innermost layers of the ITS eonsist of the Silieon Pixel Deteetor (SPD), whieh provides 
two layers of silieon pixel sensors at radii 3.9 em and 7.6 em from the beam axis and eovers the 
full azimuth over Ip | < 2 and |p | < 1.4, respeetively. The eombined information of the ITS and 
TPC ean determine the momenta of eharged partieles from low momentum (px ~ 0.15 GeV/c) 
to high momentum (px ~ 100 GeV/c) in |p | < 0.9 and full azimuth. 

The EMCal is a Pb-seintillator sampling ealorimeter, whieh eovers 107 degrees in azimuth and 
Ip I < 0.7. It eonsists of 10 supermodules with a total of 11520 individual towers eaeh eovering 
an angular region Ap x Acp = 0.014 x 0.014 whieh are read out by avalanehe photodiodes. 

The data were reeorded in 2011 for Pb-Pb eollisions at /^'nn = 2.76 TeV using a set of triggers 
based on the hit multiplieity reeorded by the VZERO deteetor, whieh eonsists of segmented 
seintillators eovering the full azimuth over 2.8 < p <5.1 (VZERO-A) and —3.7 < p < —1.7 
(VZERO-C). 

3 Data analysis 

A total of 11.5M (15/ib^/ and 5.7M (3.7/ib^/ events with VZERO multiplieities eorrespond- 
ing to 0-10% and 10-30% most eentral events were seleeted using the eentrality determination 
as deseribed in Ref. [fSOll . The aeeepted events, reeonstrueted as deseribed in Ref. lISTlI . were 
required to have a primary reeonstrueted vertex within 10 em of the eenter of the deteetor. 

Reeonstrueted traeks were required to have at least 3 hits in the ITS used in the fit to ensure ad¬ 
equate traek momentum resolution for jet reeonstruetion. For traeks without any hit in the SPD, 
the primary vertex loeation was used in addition to the TPC and ITS hits for the momentum 
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determination of the traek. This redueed the azimuthal dependenee of the traek reeonstruetion 
effieieney due to the non-uniform SPD response, without ereating traek eolleetions with dras- 
tieally differing momentum resolutions. Aeeepted traeks were required to be measured with 
0.15 < pt < 100 GeV/c in |t]| < 0.9, and to have at least 70 TPC spaee-points and no less 
than 80% of the geometrieally findable spaee-points in the TPC. The traeking effieieney was 
estimated from simulations of the deteetor response using GEANT3 Il52l with the HIJING [[5^ 
event generator as input. In 0-10% eollisions, it is about 56% at 0.15 GeV/c, about 83% at 
1.5 GeV/c and then deereases to 81% at 3 GeV/c, after whieh it inereases and levels off to about 
83% at above 6.5 GeV/c. In 10-30% eollisions, the traeking effieieney follows a similar 
dependenee pattern, with absolute values of the effieieney that are 1 to 2% higher eompared 
to 0-10% eollisions. The momentum resolution dpj/pj, estimated on a traek-by-traek basis 
using the eovarianee matrix of the traek fit, is about 1% at 1.0 GeV/c and about 3% at 50 GeV/c. 
Traeks with pj > 50 GeV/c were only a small eontribution to the inelusive jet population eon- 
sidered in this analysis, for example only 20% of the jets with pj^ jet larger than 100 GeV/c were 
found to eontain a traek above 50 GeV/c. 

EMCal eells with a ealibrated response of more than 50 MeV were elustered prior to inelusion 
in the jet finder by a elustering algorithm whieh required eaeh eluster to only have a single lo- 
eal maximum ||49l . Interaetions of slow neutrons or highly ionizing partieles in the avalanehe 
photodiodes ereate elusters with large apparent energy, but anomalously small number of eon- 
tributing eells, and are removed from the analysis. A non-linearity eorreetion, derived from 
eleetron test beam data, of about 7% at 0.5 GeV and negligible above 3 GeV, was applied to the 
elusters’ energies. The energy resolution obtained from eleetron test beam data is about 15% at 
0.5 GeV and better than 5% above 3 GeV. 

Unlike eleetrons and photons, whieh deposit their full energy in the EMCal via eleetromagnetie 
showering, eharged hadrons deposit energy in the EMCal, mostly via minimum ionization, but 
also via nuelear interaetions whieh generate hadronie showers. To avoid double eounting, the 
energy deposited in the EMCal by eharged partieles that were already reeonstrueted as traeks, 
the elusters’ energies were eorreeted by the following proeedure [|5^ : All traeks with pi > 
0.15 GeV/c were propagated to the average eluster depth within the EMCal, and then assoeiated 
to elusters with Ej > 0.15 GeV within the window |A? 7 | < 0.015 and |A(p| < 0.025. Traeks 
were always matehed to their nearest eluster, while elusters were allowed to have multiple traek 
matehes. Clusters with matehed traeks were eorreeted for eharged partiele eontamination by 
removing the fraetion / = 100% of the sum of the momenta of all matehed traeks from the 
eluster energy, as done in [[54l . Clusters with Ej > 0.30 GeV after this eorreetion were used in 
this analysis. 

The eolleetion of traeks and eorreeted EMCal elusters was then assembled into jets using the 
anti-kj or the kj algorithms in the EastJet paekage Il55ll with a resolution parameter of R = 0.2. 
Only those jets that were at least R away from the EMCal boundaries of | t] | <0.7 and 1 .4 < (^ < 
n, and thus fully eontained within the EMCal aeeeptanee, were kept in the analysis whieh limits 
the effeet of the aeeeptanee boundaries on the measured jet speetrum. Jets reeonstrueted by the 
anti-kj algorithm were used to quantify signal jets, while jets reeonstrueted by the kj algorithm 
were used to quantify the eontribution from the underlying event. 

The signal speetrum formed from the reeonstrueted jets is affeeted by the eontribution from the 
underlying event. In order to suppress the eontribution of the baekground to the measurement of 
the jet energy, we followed the approaeh deseribed in Refs. [|45ll46l . whieh addresses the aver- 
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Fig. 1 : The 5pj distribution for R = 0.2 with the random-cone and the embedded-track methods in the 
10% most central events, with = 60 GeV/c for the embedded-track method. 


age additive contribution to the jet momentum on a jet-by-jet basis. The underlying background 
momentum density was estimated event-by-event using the median of PTjetMjet’ where 
the uncorrected energy and Ajet is the area of jets reconstructed with the kj algorithm. Due to 
the limited acceptance of the EMCal, Pch, the median of the event-by-event momentum density 
distribution obtained from charged jets (i.e. jets reconstructed from tracks only) in |T]jet| < 0.5 
and full azimuthal acceptance was used. Then, Pscaied was determined by scaling p^h using 
a centrality-dependent factor. This factor is obtained from a parametrization of the measure¬ 
ment of the charged-to-neutral energy ratio, using tracks and corrected clusters in the EMCal 
acceptance. In 0-10% central Pb-Pb collisions, the average charged background momentum 
density was (Pch) ~ 110 GeV/c. After scaling to include the neutral component we obtained 
(pscaied) ~ 190 GcV/c, wMch Corresponds to an average contribution of the underlying event of 
about 24 GeV/c in a cone ofR = 0.2. In 10-30% central Pb-Pb collisions (pscaied) decreases 
to ~ 130 GeV/c. Eor every signal jet reconstructed with the anti-kj algorithm, the background 
density scaled by the area of the reconstructed signal jet was subtracted from the reconstructed 
transverse momentum of the signal jet according to Pj“t = Px^et ~ Pscaied '^Ijet- 

Region-to-region background fluctuations lead to a smearing of the reconstructed jet energy. 
Their magnitude was estimated as described in Refs. P3] |4^ in two different ways: (1) by 
taking the scalar sum of the px of all particles found in a cone randomly placed in the event, 
referred to as random-cone method, and (2) embedding a single particle in the event and inspect- 
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ing the anti-A:T jet that eontains that embedded partiele, referred to as embedded traek method. 
The first method does not rely on any assumptions about the strueture of the baekground itself 
and gives approximately the same baekground fluetuation as embedding a traek with infinite 
momentum for anti-kx jets. The seeond method should be able to reproduee the baekground 
as seen by the anti-kx algorithm more direetly. The baekground fluetuations were quantified 
by 5/)x = Px°"^ “ Pscaied ' for the random-eone method, and 8pj = Px jet ~ for the 

embedded-traek method with a minimum of = 10 GeV/c for the px of the embedded 
traek. Above 10 GeV/c the resulting 5px distribution does not depend on the px of the embed¬ 
ded partiele. The 5px distributions for the two methods in the 10% most eentral eollisions are 
shown in Fig. [I for pj°^^ = 60 GeV/c. The two methods appear to provide the same quantitative 
response to the baekground fluetuations, with only marginal differenees mainly due to small jet 
area fluetuations in the embedding traek method. The widths of the d px distributions are about 
6 GeV/c. The left-hand-side (LHS) of the distribution is Gaussian-like and is dominated by 
soft partiele produetion. To determine its width, the distributions were fitted reeursively with a 
Gaussian funetion in the range -|- using the mean and width of the 

5px distribution as starting values for o and p. The LHS width is about 5 GeV/c in 0-10% and 
about 3.5 GeV/c in 10-30% events. The right-hand side has additional eontributions from hard 
seattering proeesses, and the resulting non-Gaussian tail at high 5px is due to overlapping jets. 
The random-eone method was used as the baseline in this analysis for ereating the response 
matrix used in unfolding, while the single partiele embedding method was used to study the 
sensitivity of the results to the method. 

Additionally, signal jets were required to eontain a eharged traek with a transverse momentum 
of at least 5 GeV/c and a minimum baekground subtraeted Px jet ‘^f GeV/c for 0-10% and of 
20 GeV/c for 10-30% most eentral events, whieh roughly eorresponds to 5 (7 of the 5px dis¬ 
tribution, in order to suppress the eontribution of eombinatorial jets, i.e. from jets reeonstrueted 
mainly from upward fluetuations of the soft-partiele baekground. 

Both the average baekground and the baekground fluetuations are averaged over all possible 
orientations of the event plane, namely it is assumed that the signal jet sample being analyzed 
is isotropieally distributed with respeet to the event plane. However, the jet sample may show 
some degree of eorrelation with the event plane, both for physieal reasons (e.g. path length 
dependenee of jet energy loss) or as a result of the euts applied in the analysis (most notably the 
requirement on the leading hadron px). Sinee the baekground is also eorrelated with the event 
plane due to flow (v 2 ) ffTOl . a question may arise about the validity of this approaeh. Upper 
limits on the magnitude of these effeets have been estimated by using random eones biased 
towards the event plane, either by requiring the presenee of a 5 GeV/c traek or by weighting the 
distribution using an upper limit on the jet vx of 0.1. In both oases, the upper limits on the shift 
of the jet energy seale (JES) were found to be smaller than 0.1 GeV/c. 

4 Unfolding 

The measured jet speetra are distorted by the response of the deteetors used in the measurement 
and the baekground fluetuations in the underlying event. To oorreot for these effeets we used 
an “unfolding” proeedure, as desoribed in Ref. Il46ll . The oorreoted distribution p^J^jgj and the 
measured distribution Pj “j are related by a eonvolution through the response matrix RMtot = 
RMbkg X RMdet, where RM^et parametrizes the deteetor response and RM^kg the baekground 
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fluetuations. The unfolding proeedure operates under the assumption that = RMtot x Pj“gf. 
Both baekground fluetuations and the deteetor response to jets are uniform within the rj and cp 
aeeeptanees, whieh is a preeondition for the faetorized approaeh used in building RMtot- 

The deteetor response for jet reeonstruetion was obtained using pp events simulated with the 
PYTHIA6 ll5^ event generator (tune A [|57ll '). Jets were reeonstrueted both at “generator level” 
and at “deteetor level” using the anti-kj algorithm. Generator-level simulations utilized only 
prompt partieles originating from the eollision (with cT < 1 em), direetly from the event gen¬ 
erator output, without aeeounting for deteetor effeets; deteetor-level simulations also ineluded 
a detailed partiele transport and deteetor response simulation based on GEANT3 liS^ with the 
deteetor response set to the Pb-Pb eonfiguration. During deteetor-level jet reeonstruetion, an 
additional px-dependent traeking ineffleieney was introdueed in order to aeeount for the larger 
ineffieieney due to the larger oeeupaney effeets in eentral Pb-Pb events eompared to pp events. 
Oeeupaney effeets have been estimated eomparing the traeking performanee in PYTHIA and 
HIJING simulations, whieh represent pp and Pb-Pb events [l53ll . The oeeupaney effeets in 
eentral HIJING events are larger for px < 0.5 GeV/c where the effleieney is about 4% lower 
eompared to PYTHIA, and then levels off to about 2% lower for px > 2 GeV/c. In semi-eentral 
HIJING events, oeeupaney effeets on the traeking effleieney amount to no more than 2% at low 
Px and about 1% for px > 2 GeV/c. Other than this traeking effleieney eorreetion, the detee¬ 
tor response to jets was assumed to be the same in Pb-Pb events as in the PYTHIA simulated 
pp eollisions. 

The generator-level and deteetor-level jets were matehed based on the Euelidean distanee be¬ 
tween their jet axes in pseudorapidity and azimuthal angle. It was ensured that the matehing 
operation is bijeetive: eaeh generator-level jet was matehed to at most one deteetor-level jet fl6l . 
Every matehed jet pair eorresponds to an entry in the deteetor response matrix, RM^et- An un- 
matehed generator-level jet represents a jet that was not reeonstrueted, and this distribution was 
used to determine the jet reeonstruetion effleieney. In 0-10% Pb-Pb events, the deteetor jet 
reeonstruetion effleieney was found to be 90% at 40 GeV/c and 95% above 70 GeV/c, lim¬ 
ited mainly by the traek reeonstruetion effleieney of the leading eharged partiele. As deseribed 
above, at deteetor level the eonstituent eut was 150 MeV/c for traeks, and 300 MeV for elusters 
after the eluster energy is eorreeted for eharged partiele energy eontamination. However, at 
generator level no sueh out is applied, and henee the reeonstrueted jets are eorreeted to a oon- 
stituent eharged partiele momentum of 0 MeV/c and to a eonstituent eluster energy of 0 MeV 
in the unfolding prooess. A net negative shift of the JES at deteetor level was obtained, whieh 
originates mainly from traeking ineffleieney and unreoonstruoted partieles, sueh as neutrons and 
though the subtraotion proeedure for energy deposits by eharged partieles in the EMCal and 
missing seoondary partieles from weak deoays eontribute to the shift [f54l . The JES eorreetion 
applied through the response matrix is about 23% at Px’jet GeV/c and 29% at 120 GeV/c 
independent of eentrality. 

The RMbkg matrix was eonstrueted row-by-row by taking the dpj distribution and shifting it 
along the Px jet amount Px'jet corresponding to eaeh row (Toeplitz matrix). This 

matrix eonstruetion method assumes that the response of the jet speetrum to baekground flue- 
tuations is independent of the jet momentum. 

The px-dependenee of the jet momentum resolution 0'(pj “J /Px^gj is different for the baek¬ 
ground and deteetor eontributions ll4i^ . The eontribution from baekground fluetuations is dom¬ 
inant at low Px’jet proportional to 1/p^gj, whereas the eontribution from deteetor ef- 
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feets is fairly eonstant with The eross-over between the two eontributions happens at 

pJPjgt ~ 30 GeV/c. The eombined jet momentum resolution is about 23% at Pj“gj of 40 GeV/c 
and 20% at 120 GeV/c for 0-10% eollisions, while it is 24% at p^^gj of 30 GeV/c and 20% at 
100 GeV/c for 10-30%. 

Two unfolding algorithms with different regularization proeedures were used for eorreeting the 
measured jet speetrum: the 'y} minimization method [|58ll with a log-log-regularization and 
the generalized Singular Value Deeomposition (SVD) method [f59ll . as implemented in RooUn- 
fold [[60l . whieh was used for the default value of the data points. The measured speetrum 
used as an input to the unfolding was in the range 30 < px, jet < 120 GeV/c for 0-10% and 
20 < px, jet < 100 GeV/c for 10-30% eollisions. A smoothed version of the measured speetrum 
was used as the prior, so that the statistieal fluetuations within the data were not magnified in the 
unfolding proeess. The regularization parameter used for SVD unfolding is A: = 5. The value 
of k is ehosen sueh that it eorresponds to the d veetor magnitude of 1, and Pearson eoeffieients 
whieh do not show a large variation in the eorrelation between neighboring j!?x bins. 

The eorreeted jet speetra are reported for 40 < Px, jet < 120 GeV/c in 0-10%, and for 30 < 
Pt, jet < 100 GeV/c in 10-30% where the effieieney due to these kinematie euts is high, ap¬ 
proximately 90%. It was verified that the eut on the reeonstrueted jet px has a negligible effeet 
in the reported px region of the final result, as long as the requirement on the leading eharged 
traek px is at least 5 GeV/c. If this threshold is redueed, the eut on the minimum reeonstrueted 
jet Px beeomes erueial for unfolding stability. 

The analysis proeedures in the 10% most eentral eollisions were tested with two different Monte 
Carlo (MC) models, where events were eonstrueted by embedding jets into a soft baekground. 
The first test verified the robustness of the unfolding framework with the inelusion of fake 
“jets” that are elustered from the soft baekground, whieh did not originate from a hard proeess. 
The seeond model tested the assumption that the baekground and deteetor responses ean be 
faetorized. 

In the first model, the soft baekground of both eharged and neutral partieles was modeled with 
3100 < A^tracks < 5150 where the partiele transverse momenta were taken from a Boltzmann 
distribution with a temperature of 550 MeV. This model ereated a fluetuating baekground sim¬ 
ilar to that of the 0-10% Pb-Pb data; e.g. the baekground fluetuations, as estimated via the 
5px distributions, eoineide within few pereent. Jets were reeonstrueted at generator level in 
PYTHIA-only events and at deteetor level, with the added baekground. The first model val¬ 
idated the baekground subtraetion teehnique, and in partieular the stability of the unfolding 
method against the eontribution from the residual eombinatorial baekground. In the seeond 
model, the baekground was taken from real 0-10% Pb-Pb events. The eharged partiele eor- 
reetion for the EMCal elusters was applied after embedding. Only jets with at least 1 GeV/c 
of transverse momentum eoming from the embedded PYTHIA event were seleeted for the test. 
This is needed to rejeet the signal from hard seatterings in the data, but also removes most of 
the eombinatorial jets from the Pb-Pb underlying event. The seeond model was used to test 
the validity of the eharged partiele eorreetion applied to the EMCal elusters, in partieular in 
the interplay between the underlying event and the jets. It also validates eertain aspeets of the 
eorreetions applied for the baekground fluetuations, e.g. the unsmearing of the jet px due to 
baekground fluetuations or the overlap with low momentum jets. Baekground traeks and elus¬ 
ters eould be matehed to jet traeks and elusters or viee versa, so that the eorreetion for eharged 
partiele eontamination eould potentially eause an over-subtraetion that is not eorreeted for in 
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the unfolding proeedure. These Monte Carlo tests showed that the analysis proeedures outlined 
above, ineluding unfolding, reeovered the input speetrum within the statistieal and systematie 
uneertainties of the models. 


Category Relative uneertainty (%) 



„min 
Pt, jet 

y,max 
^T, jet 

Min. 

Max. 

Avg. 

Tracking efficiency 

7.7 

11.3 

7.3 

11.3 

8.8 

Track momentum resolution 

1.0 

1.0 

1.0 

1.0 

1.0 

Charged particle correction 

0.7 

2.7 

0.7 

6.4 

3.7 

EMCal clusterizer 

3.2 

1.8 

0.1 

3.2 

1.4 

EMCal response 

4.4 

4.4 

4.4 

4.4 

4.4 

Background fluctuations 

3.9 

2.7 

2.3 

3.9 

2.8 

Jet raw pj cuts 

2.6 

6.7 

1.5 

6.7 

3.6 

Combinatorial jets 

0.3 

0.5 

0.0 

0.5 

0.2 

Total correlated uncertainty 

10.6 

14.5 

10.6 

14.5 

12.0 

Unfolding method 

0.1 

10.0 

0.1 

15.5 

6.6 

SVD reg. param. k = A 

3.6 

11.7 

2.4 

11.7 

6.0 

SVD reg. param. k = 6 

7.2 

2.7 

1.5 

8.8 

5.3 

Prior choice 1 

1.9 

4.0 

0.2 

4.0 

1.6 

Prior choice 2 

2.1 

1.4 

0.1 

2.1 

0.9 

Total shape uncertainty 

3.8 

7.2 

2.7 

7.4 

5.3 


Table 1: Summary of systematic uncertainties for 0-10% most central collisions. The first column is the 
uncertainty at the minimum of 40 GeV/c, the second column is the uncertainty at the maximum 
Pr^jet GeV/c. The minimum and maximum columns give the extreme, and the last column 

gives the average systematic uncertainty over the entire pj range. The total correlated uncertainty was 
calculated by adding the components in quadrature, while the shape uncertainty was calculated as the a 
of the different variations (see text for details). 


5 Results 

The unfolded jet spectra in 0—10% and 10—30% central collisions are displayed in Fig. To 
compare the spectra with the spectrum measured in pp collisions, the yield is divided by the 
number of binary collisions, which is Acoii = 1501 ± 167 for 0-10% and 743 ±79 for 10-30% 
collisions, as estimated from a Glauber MC calculation IISOll . 

The systematic uncertainties on the jet spectrum are summarized in Tab. for the 0-10% cen¬ 
trality class. For the 10-30% centrality class the corresponding uncertainties differ, on average, 
by 2% or less. The systematic uncertainties were divided into two categories: correlated uncer¬ 
tainties and shape uncertainties. The correlated uncertainties result dominantly from uncertain¬ 
ties on the JES, such as the uncertainty of the tracking efficiency, that will shift the entire jet 
spectrum in one direction, whereas the shape uncertainties are related to the unfolding and can 
distort the slope of the spectrum. 

The dominant correlated uncertainty on the jet spectrum of about 9% arises from the uncertainty 
on the tracking efficiency. It is estimated by varying the tracking efficiency by 5% in determin¬ 
ing RM(jet and unfolding the spectrum. The uncertainty due to the correction procedure for the 
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Fig. 2 : The spectra of R = 0.2 jets with a leading track requirement of 5 GeV/c in 0-10% and 10-30% 
most central Pb-Pb collisions scaled by 1 /Ncoii and in inelastic pp collisions at .^/^nn = 2.76 TeV. The 
uncertainties on the normalization are about 11 % for the Pb-Pb data from the uncertainty on Acoii and 
about 8% for the pp data from the total inelastic cross section. 


charged particle double counting in the EMCal of about 4% was determined by varying / from 
100% to 30% in both the measured spectrum and the RM^et- The determination of the uncer¬ 
tainties from other EMCal response related uncertainties as EMCal energy scale, EMCal energy 
resolution, and EMCal non-linearity is outlined in [[541 and combined leads to an uncertainty of 
4.4%. The uncertainty arising from the choice of the EMCal clustering algorithm is determined 
by using a different clusterizing method, that forms fixed-size clusters from 3x3 towers. Eor 
the background fluctuations, the response matrix RMbkg was constructed with the single-track 
embedding method for determining dpr, as discussed above. To estimate the sensitivity of the 
unfolding to the raw jet pj selection, the pj range of input spectra is varied by extending the 
range at both the low and high ends by ±5 GeV/c. The influence of combinatorial jets, esti¬ 
mated by varying the low edge of the unfolded spectrum from 0 to up to 10 GeV/c was found 
to be negligible. Since all sources of uncertainty are independent, each contribution is added 
in quadrature to obtain the final correlated uncertainty of 10.6% to 14.5% as listed in Tab. 
The uncertainty on the JES is 2.4% to 3.2% and can be obtained by dividing the uncertainties 
listed in Tab. [^by 4.5, where the exponent n = 4.5 was obtained by fitting a power law to the 
measured spectrum. 

The shape uncertainty is dominated by the regularization used in the unfolding and can be 
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Pt,, (GeV/c) 

Fig. 3 : Ratio of the jet spectrum with a leading track pt > 5 GeV/c over the inclusive jet spectrum for 
R = 0.2 in pp collisions at ^/s = 2.76 TeV. 


divided into two eomponents: the method by whieh the solution is regularized, e.g. instead 
of the SVD unfolding, and the variation of the regularization proeess within a given method. 
The regularization is done by adding a penalty term in the method and by ignoring the 
eomponents of the SVD deeomposition that are dominated by statistieal fluetuations. For the 
SVD method, the regularization k faetor is an integer value and thus ean only be varied in integer 
steps. The uneertainty related to the ehoiee of the prior is estimated by varying the exponent 
of the power law funetion extraeted from the reeonstrueted speetrum by ±0.5, whieh is used to 
eonstruet the prior. The uneertainty related to the ehoiee of the prior is estimated by varying 
the exponent n = 4.5 by ±0.5 to seale the prior. The differenees in the unfolded speetrum with 
these variations are summarized in Tab. These variations in the regularization strategy are 
eombined assuming that they eonstitute independent measurements. The final shape uneertainty 
is thus obtained by summing them in quadrature and dividing by the square root of the number 
of variations. 

The jet speetrum in pp eollisions was measured in the same way as reported in Ref. Il54ll . but 
with the 5 GeV/c leading eharged partiele requirement neeessary for the Pb-Pb analysis. The 
resulting speetrum normalized per inelastie pp eollision is shown in Fig.|^ In order to determine 
the effeet of the leading traek requirement in pp eollisions, the ratio of the jet speetra with a 
5 GeV/c leading traek requirement (the biased jet sample), over the speetrum of jets without a 
leading traek requirement (the inelusive jet sample) with resolution parameter R = 0.2 is shown 
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in Fig.[^ Systematic uncertainties in the ratio were evaluated by removing the uncertainties that 
are correlated between the spectra obtained with and without the cut on the leading particle. As 
can be seen in Fig.j^ for pj jet above 50 GeV/c more than 95% of all reconstructed jets have at 
least one track with a pj greater than 5 GeV/c. PYTHIA tune A (but also other common tunes 
like the Perugia tunes [l57l ') accurately describe the measured ratio. 


CO 
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Fig. 4: Ratios of jet spectra with different leading track pj requirements (“0 over 5”, “3 over 5”, “7 over 
5”and “10 over 5”) for/? = 0.2 jets in 0-10% Pb-Pb collisions at = 2.76 TeV. The solid black lines 
represent predictions from PYTHIA. 


The influence of the leading track requirement in the Pb-Pb measurement, nominally set to 
5 GeV/c was tested by varying it by 40%, i.e. reducing it to 3 and increasing it to 7 GeV/c, and 
with the more extreme values of 0 and 10 GeV/c. The ratios of jet spectra with the different 
leading track pj biases, after all corrections, are shown in Fig. |^for /? = 0.2 jets in 0-10% 
central Pb-Pb collisions at = 2.76 TeV. The corrections to these different jet spectra were 
done using the same unfolding procedure as the nominal spectrum with leading track pj bias of 
5 GeV/c, with a slightly modified response matrix which accounts for the different biases. Since 
the unfolding procedure weakens the correlation between the statistical fluctuations of the jet 
spectra with different leading track requirements, the statistical uncertainties have been added 
in quadrature in the ratio. The systematic shape uncertainty is due to the unfolding procedure, 
and has been treated as completely uncorrelated in the ratio. The correlated uncertainty is 
primarily due to the uncertainty on the JES, which is highly correlated between the various 
spectra. The systematic variations in the unfolding procedure have been applied consistently 
for both the denominator (with a leading track pj > 5 GeV/c) and the numerators (with a 0, 3, 
7 and 10 GeV/c leading track bias), and the resulting difference in the ratios has been taken as 
a systematic uncertainty. The jet spectra with leading track requirements of 3 and 0 GeV/c are 
consistent with the baseline measurement with a 5 GeV/c requirement. The unfolding is not as 
stable as with a 5 GeV/c requirement, which leads to a larger systematic uncertainty due to the 
unfolding correction procedure, especially for the inclusive spectrum. All measurements of the 
ratio of jet spectra with different leading track biases, particularly those with a higher leading 
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track pj requirement than the nominal, are well described by PYTHIA 6 (tune A), within one 
sigma of the uncertainties or less. 



Fig. 5: /?aa for R = 0.2 jets with the leading traek requirement of 5 GeV/c in 0-10% (left) and 10-30% 
(right) most central Pb-Pb collisions compared to calculations from YaJEM llbTI and JEWEL |[62l . The 
boxes at 7 ?aa = 1 represent the systematic uncertainty on Taa- 


The nuclear modification factor, Raa, is defined as the ratio of the jet spectrum in Pb-Pb di¬ 
vided by the spectrum in pp collisions scaled by Nco\\. It is constructed such that Raa equals 
unity if there is no net nuclear modification of the spectrum in Pb-Pb collisions as compared 
to an incoherent superposition of independent pp collisions. The resulting Raa of jets with a 
5 GeV/c leading track requirement for R = 0.2 in the 0—10% and 10—30% central Pb-Pb col¬ 
lisions is reported in Fig. The systematic and statistical uncertainties from the Pb-Pb and pp 
measurements (see Fig.|^ are added in quadrature. The resulting uncertainty on the normaliza¬ 
tion is from scaling the pp cross section with the nuclear overlap Taa = 23.5 ±0.87 mb^^ for 
0-10% and 1 L6±0.60 mb^^ for 10-30% collisions. As can be seen, jets in the measured jet 
range are strongly suppressed. The average Raa in both 0-10% and 10-30% central events was 
found to have a negligible jet dependence. In the 10% most central events, combining the 
statistical and systematic uncertainty in quadrature, the average Raa is found to be 0.28 ±0.04. 
The suppression is smaller in magnitude in the 10-30% central events, leading to an average 
J?AAof 0.35 ±0.04. These results qualitatively agree with the suppression obtained from mea¬ 
surements using charged-particle jets [|46]l . though the jet energy scale is not the same in both 
cases, and so a direct comparison is not possible. Furthermore, the results are consistent with 
the Raa reported by ATFAS for R = 0.4 jets scaled by the ratio of the yields with the different 
resolution parameters in different px, jet bins [|3^I4T]| . 

In order to interpret the results and move to a more quantitative understanding of jet quenching 
mechanisms, a comparison of the measured Raa in 0-10% central collisions to calculations 
from two different models is also shown in Fig. The first model, YaJEM [IMI . uses a 2-i-lD 
hydrodynamical calculation and a Glauber MC for the initial geometry, as well as a FO pQCD 
calculation to determine the outgoing partons. Parton showers are modified by a medium- 
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indueed inerease of the virtuality during their evolution through the medium. The Lund model in 
PYTHIA is used for hadronization into final state particles. The kinematics of the virtual partons 
in the evolving partonic shower were modified with a parameter related to the two transport 
coefficients, q and e, that describes how strongly a parton of a given momentum couples to the 
medium. The parameter was fixed so that the model accurately describes the 7 ?aa for charged 
hadrons at 10 GeV/c IfTTlI . but no additional changes were made for the prediction of the jet 
7 ?aa- The second model, JEWEL [[62||, takes a different approach in the description of the 
parton-medium interaction by giving a microscopical description of the transport coefficient, 
q. Essentially each scattering of the initial parton with medium partons is computed and the 
average over all scatters determines q. JEWEL uses a combination of Glauber and PYTHIA to 
determine the initial geometry, a ID Bjorken expansion for the medium evolution, and PYTHIA 
for hadronization into final state particles. The transverse medium density profile in JEWEL is 
proportional to the density of wounded nucleons combined with a ID Bjorken expansion for the 
time evolution. Hard scatters are generated according to Glauber collision geometry, and suffer 
from elastic and radiative energy loss in the medium, including a Monte Carlo implementation 
of LPM interference effects. PYTHIA is used for the hadronisation of final state particles. 
Despite their different approaches, both calculations are found to reproduce the jet suppression. 
YaJEM, however, exhibits a slightly steeper increase with jet pi than the data. The calculated 
are 1.690 for YaJEM and 0.368 for JEWEL, obtained by comparing the models with the 
data. Additional measurements will be needed in order to further constrain the models, such 
as measuring the jet suppression relative to the event plane angle, which would require a more 
accurate modeling of the path-length dependence of jet quenching. 

6 Summary 

The transverse momentum (px) spectrum and nuclear modification factor (Raa) of jets recon¬ 
structed from charged particles measured by the ALICE tracking system and neutral energy 
measured by the ALICE Electromagnetic Calorimeter are measured with R = 0.2 in the range 
of 40 < px, jet < 120 GeV/c for 0-10% and in 30 < px, jet < 100 GeV/c for 10-30% most cen¬ 
tral Pb-Pb collisions at = 2.76 TeV were measured. The jets were required to contain at 
least one charged particle with px > 5 GeV/c. The effect of this requirement on the reported 
Raa was evaluated by the ratios of the jet spectra with the 5 GeV/c to no requirement compared 
to expectations on PYTHIA, and found not to have an observable effect within the uncertainties 
of the measurement. Jets with 40 < px, jet < 120 GeV/c are strongly suppressed in the 10% 
most central events, with Raa about 0.28 ± 0.04, independent of px. jet within the uncertainties 
of the measurement. The suppression in 10-30% events is 0.35 ±0.04, slightly less than in the 
most central events. The observed suppression is in fair agreement with expectations from two 
jet quenching model calculations. 
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